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**F ILE**1D**OTSPOWCGJ 

000000 TITTTTITTT SSSSSSSS PPPPPPPP 000000 wie wid cccccccc GGGGGGGG JJ 
000000 TITTTTTITT SSSSSSSS_ PPPPPPPP 000000 wd Wel ceccccce GGGGGGGG JJ 

00 00 TT $$ PP PP 00 00 WW ww CC GG JJ 

00 00 TT SS PP PP 00 00 WwW ww CC GG JJ 

00 00 TT SS PP PP 00 00 WW ww CC GG JJ 

00 00 TT Ss PP PP 00 00 WW ww CC GG JJ 

00 00 TT SSSSSS PPPPPPPP 00 00 WW ww CC GG JJ 

00 00 TT SSSSSS PPPPPPPP 00 00 WW ww CC GG JJ 

00 00 TT SS PP 00 00 WwW Ww Ww CC GG GGGGGG JJ JJ 

00 00 TT SS PP 00 00 WW Ww Ww CC GG GGGGGG JJ JJ 

00 00 TT SS PP 00 00 WwWW WhWWw CC GG GG JJ JJ eoee 

00 00 TT SS PP 00 WhWW Wwww CC GG GG JJ JJ eee 
000000 TT SSSSSSSS PP 000000 ‘ells wid ccccccce GGGGGG JJJJJJ cece 
000000 TT SSSSSSSS PP 000000 wis we ccccccce GGGGGG JJJJdJ eeee 

LL III SSSSSSSS 

LL HI SSSSSSSS 

LL I] SS 

LL I] $$ 

LL I] SS 

LL II SS 

LL I] SSSSSS 

LL I] SSSSSS 

LL I] SS 

LL I] SS 

LL I] SS 

LL I] SS 

LELLLLLLLL III] SSSSSSSS 

LLELLLLLLLL III! SSSSSSSS 
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MEUM OOD NAULS WN (SO ODONAUES WN $9 OONAUES WN) OWONOUS wn 
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0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
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COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


v ® 
® ® 
® ® 
& ® 
a ® 
* THIS SOFTWARE 1S FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
* ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
* INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER * 
* COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ARY 
* OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 
:® TRANSFERRED. * 
® a 
* THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE * 
* AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT + 
* CORPORATION. ® 
& * 
*® ® 
** ® 
*® * 
*® ® 
‘ ® 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 


PRRRRARAAALAAAALLAAA ALAS ALA SE ALES ASE LESSEE EE ESE RSE EERE RSE REE REE REET R EE SY 


FACILITY: Language support library - user callable 
“ABSTRACT: 


G COMPLEX*16 base to arecEnss power. 
Floating overflow can occu 

Undefined Sgpgnent tet ton can occur if 
base = (0.,0.) and exp <=0 


VERSION: 1 


HISTORY: 
AUTHOR: 
Steven B. Lionel, 27-July-1979 


TSSPOWCGJ - G COMPLEX*16 ** INTEGER*4 o a. atin 16-SEP-1984 01:56:33 VAX/VMS Macro V04-00 P 
a HISTORY “© Detailed Current Bate History S-SEPo19RG $1:39;8% LMTAKTE. Sarsorseouccy mars1 2% 3 
2 ~SBTTL HISTORY ; Detailed Current Edit History 
i 
; dit History 
: mt - Adapted te OTS$POWCJ vereten 1-003. SBL 27*July-1979 
: 1-002 - a6 x bug 3 965°8ks Ng Mert ined exponentiation with negative powers. 
: 1-003 - ise general mode slircotian. SBL 30-Nov-1981 
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TSSPOWCGJ - G COMPLEX*16 ** INTEGER*®4 power routin 16-SEP-19 1:56:33 VAX/VMS Macro v04-00 
a ts DECLARATIONS . aaa 0:39:33 MTHRTL. ga ¢40 OTSPOWCGJ.MAR;1 
He ‘ -SBTTL DECLARATIONS 
8 § ; INCLUDE FILES: 
62 ; EXTERNAL SYMBOLS: 
38 8§ 3 
0 64 
00 65 -DSABL GBL 
00 66 EXTRN MTHSSSIGNAL 3; Math error routine 
4 67 -EXTRN OTSSDIVCG_R3 ; COMPLEX division routine 
000 68 eEXTRN MTHS$K_UNDEXP 
0000 $? 
0000 0; 
0000 71 =; MACROS: 
0000 if $ 
0000 7 
0000 74 ; 
0000 75 ; EQUATED SYMBOLS: 
0000 76 ; 
0000 77 
0000 78 ; 
0000 79 ; OWN STORAGE: 
0000 80 ; 
0000 81 
0000 Hf 3 
0000 83 ; PSECT DECLARATIONS: 
0000 84; 
0000 85 
00000000 86 -PSECT _OTSSCODE PIC,SHR,LONG, EXE, NOWRT 
464 rH 3; program section for OTS$ code 
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T 


C1 | 
OTSSPOWCGJ - G COMPLEX*16 ** INTEGER*4 power routin 16-SEP-1984 01:56:33 VAX/VMS Macro v04-00 e 4 
1-003 OTSSPOUCED RS = G COMPLER STO oe INTERES CoSEP TORS $1258 Tee EAU eae Oe 0 many 2% cd 
0 LSBTTL OTSSPOWCGJ_R3 = G COMPLEX®16 ** INTEGER®4 


poe 
; FUNCTIONAL DESCRIPTION: 


G COMPLEX*16 result = G COMPLEX®*16 base ** INTEGER*4 exponent 
The COMPLEX result is given by: 


0000 
Wwrn—o 


base exponent result 

any >0 PRODUCT (base * 2**i) where 
i is each non-zero bit in 
exponent. 

(0., 0.) <=0 Undefined exponentiation. 

not (0., 0.) <0 PRODUCT (base * 2**i) where 
i is each non-zero bit in 
‘exponent:. 

not (0., 0.) =0 (1.0, 0.0) 


] 
| 
| 

Floating overflow can occur. 
Undefined exponentiation occurs if base is 0 and 
exponent is 0 or negative. 

CALLING SEQUENCE: 

result.wgc.v = OTS$POWCGJ_R3 (base.rgc.v, exponent.rl.v) 
INPUT PARAMETERS: 

base = 4 ; G COMPLEX*16 base passed by VALUE ! 

exponent = 20 ; Longword integer exponent by 


IMPLICIT INPUTS: 
NONE 


00000004 


00000014 value. 


SOoOSCSSCSOSOSOOSCSOSOSOSOSOSOSOSOOSOSOOOSOSOOOOSOoOO 
SOoOSSSOSCSOSOSOSOSOOSOOCOOSOSOOOSOOOSCOOSoOSoOO 


OUTPUT PARAMETERS: 
NONE 

IMPLICIT OUTPUTS: 
NONE 


FUNCTION VALUE: 


THE G COMPLEX*16 result is returned in registers RO-R3. 
This is a violation of the VAX calling standard, but is 
excused for compiled code support routines. 


SIDE EFFECTS: 


Modifies registers RO-R3! 

SS$_FLTOVF = A gver ites ’ 
SIGRALS MTH$_UNDEXP (82 = * UNDEFINED EXPONENTATION’) if 
base is 0 and exponent is 0 or negative. 


Sssssssssssssssssssssssssssssssssssssssssssssssssssssssss 
SOOCOOCOCOCOSOOOCOOCOOSOSOSOSSOSSSOSSSOSOCOSOCOSOSCOOSOSOOCOOOOOOOOOOOOOOOOOoOO 


BIE ARI $9 ODNAUE WIN $$ O OD NAUE WN 9 OD NAME WN 0 ODNOAU SWI OOOnNO 


a a a ak ak et ot — 8 = 8 = 9 = = = 4 = — = = —  — —) — —) — sd —)  —s 


SOOO OCOCOSOOSOSCOOOOoODOoOSoSoooO 


SOooooooooooooooooooo 


WADYVIMOW1O YS PP 
7 
ww 


WO 08 SWIM roo OONASLOWOOM 


MUIMIVITICOUNIO 


UNDEF INED: 
ASHQ 


RININPONINS @ PPP? OM 9 Ps OP Ss 9 sn in nn in 8 a tn at ts ss ss — — — — — ss —s ss 2s -s 2 2 5 2 2 5 

DBDOOOO OOOO O0 OOOO O09 6960 00 09 09 09 09 09 09 SINISE IPA AAA AAAA A UMIUIUIITUTUTUT 

FANS DOONAN WIN O ODO NAUE WIN $$ O OO NA UENO OD NAME WN OOONAUE WH" O00 
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sENTRY OTSSPOWCGJ_R3, “M<R4,R5,R6,R7,R8B> 

; disable integer overflow 
MOVa base(AP), R4 ; R4-R7 gets COMPLEX base 
MOVO base+B(AP), R6 
MOVL exponent (AP), R8 


Re = longword exponent 


GEQ 3 = | exponent : 
MNEGL Be. Re 
BBCC #0, R8, EVEN $ 
MOVG R4, R : = initial result 
MOVG R6, R 


R8 = unsigned_exponent / 2 
done if exponent was 1 
else use rest of exponent 


ROTL #1, RB, RB 
DONE 
BRB SQUAR1 


@o 
m 
© 
fo 
. 


MOVG #1, RO 


ROTL #-1, R8, RB 
BNEQ pune 


RO-R3 = initial result 
(1.0, 0.0) 


R8 = unsigned_exponent /2 
branch if exponent not 0 


Bete Se Ge Ge Ge Ge te 


] 
| 
proneh if even and clear low bit | 
RO=R 


TSTG exponent was 0, text RP(base) 
BNEQ DONE done if non-0, answer is 1.0 
TSTG R6 IP(base) better not be zero 
BNEQ DONE it isn't return 1.0 

#15, #1, RO 3; return RO-R3 = reserved operands 


ASHQ #15, » R2 
MOVZBL #MTHSK_UNDEXP, -(SP) FORTRAN error number 
CALLS #1, G*ATHSSSIGNAL vert to 32-bit condition code 


; con 
and SIGNAL MTH$_UNDEXP 


RET 


ASHL = #1, RB, RB ; R8 = ireduced exponent: / 2 
= square current base 


MULG3 R4, R6, -(SP) : (SP) = tmp = RP(base)*IP(base) 
MULG R4, RG 3; R4-R5 = RP( base) **2 

MULG R6, R6 3: R6-R7 = me ge 

SUBG R6, R4 3: R4-R5S = RP(base)*#2 - IP(base)**2 
ADDGS (SP), (SP)+, R6 ; R6-R7 = 2*(RP(base) *IP (base) ) 
BLBC RB, SQUAR ; branch if next exponent bit is 


= partial result * current power of base 


MULG3 RO, R6, -(SP) : (SP) mp = RP(part) * IP(base) 

MULG 7 R ; RO-R1 = RP(part) * RP(base) 

MULG R2, R6, -(SP) ; (SP) = tap = IP(part) * IP(base) 

SUBG2 (SP)+, RO + RO-R1 = RP(part) *RP(base)-IP(part)*I1P(base | 
MULG RG 3 R2-R3 = IP(part)*RP( base) 

ADDG2 (SP)+, R2 + RO-R3 = IP(part) *RP(base) +RP(part) *1P(base 
ASHL #-1, RB, RB ; RB = treduced exponent: / 2 


BNEG § SQUAR1 
TSTL exponent (AP) 


loop if more exponent bits left 
test exponent sign 


— 
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* ee ; done if positive 
R 3; test RP(result) 
RECIP : if non-0, OK to take reciprocal 
; RP(result) was 0, test IP(result) 
UNDEF INED 3 undefined (0.0+40.0i1) ** -n 
Ree -(SP) 3; second arg pair is divisor 
R (SP) 
-(§P) ; push (1.0,0.0) on stack 
#1, -(SP) 
#8, G*OTSSDIVCG_R3 ; RO-R3 = reciprocal 


; result in RO-R3 


Symbol table 


BASE = 

EVEN fr 8 
EVEN ¢ R 
EXPONENT = 4 
MTHSSSIGNAL eececeee § Xx 

AT UNDE eeereere x 


H x 
OTSS$OTVCG a. aereeere § =X 


St tI OOO 


RECIP 
SQUAR 
SQUAR1 
UNDEF INED 

peses see eee ee oceans + 

! Psect synopsis ! 

bo wana coenmneonen eo} 
PSECT name Allocation PSECT No. Attributes 
SS ,. 00000000 <( 0.) OO ¢ QO.) NOPIC USR CON ABS L 
-OTSSCODE 00000085 (¢ 181.) 01 ¢ 1.) PIC USR CON REL LCL 

gece ewer ese omc emer ener eo c= ot 
H Performance indicators ; 
Phase Page faults CPU Time Elapsed Time 
Initialization 30 00:00:00.10 00:00:01.0 
Command processing 108 0:00: 3-38 0:00:02.7 
Symbol tabl t % 6; : af 0:00;0 bi 
e sor Hy : . : : . 

Pass 2 5 0:00:00:51 09:00:0 ay 
Symbol table output 8 :00:00.01 0:00:00.41 
Psect synopsis output 0:00:00.03 00:00:00.06 
Cross-reference output es Sy OF 8808-80 
Assembler run totals 271 0:00:01.81 0:00:09.7 


[he working get Limit was 900 pages. 
3177 bytes (7 pages) of virtual memory were used to buffer the intermediate code. 
There were 10 pages of symbol table space allocated to hold 13 non-local and 1 local symbols. 
source Lines were read in Pass 1, producing 11 object records in Pass 2. 
pages of virtual memory were used to define 0 macros. 


+ + 
H Macro library statistics : 


Macro Library name Macros defined 
~$255$DUA28: CSYSLIBISTARLET.MLB; 2 0 
0 GETS were required to define 0 macros. 


There were no errors, warnings or information messages. 
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CL wee = NORD NOWRT NOVEC BYTE 
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s EXE RD NOWRT NOVEC LONG 


OT! 
Tal 
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GLENTASEro mn statistics” * COPLEIOG e+ INTEGERS pover'rovein AEE IRES FL:3E WARE SRS ARISE ass "OP 


MACRO/ENABLE=SUPPRESS1ION/D1 SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:OTSPOWCGJ/OBJ=0BJ$:OTSPOWCGJ MSRC$:OTSPOWCGJ/UPDATE=(ENHS$:OTSPOWCGJ) 


244 AH BT13A-SE : ENT CORPORATION 
VAX/VMS V4.0 ND PROPRIETARY 


